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Equations  and  curves  are  presented  for  estimating  height  and  site  index  of  Douglas-fir 
(Pseudotsuga  menziesii  (Mirb.)  Franco)  on  hot,  dry  sites  in  the  Willamette  National 
Forest  in  western  Oregon.  The  equations  are  based  on  the  dissected  stems  of  27  trees. 
The  curves  differ  from  those  previously  published  for  Douglas-fir.  Instructions  are 
presented  for  their  use. 

Keywords:  Increment  (height),  site  index,  Douglas-fir,  Oregon  (Willamette  National 
Forest). 

Research  Summary      This  study  provides  height  growth  and  site  index  curves  for  Douglas-fir  (Pseudotsuga 

menziesii  (Mirb.)  Franco)  on  hot,  dry  sites  in  the  Willamette  National  Forest  in  western 
Oregon.  Stems  of  40  trees  were  dissected;  27  of  the  trees  were  suitable  for  construction 
of  height  growth  and  site  index  curves,  and  they  provided  505  observations  of  height, 
site,  and  age. 

Smoothed  height  growth  curves  developed  from  these  data  were  compared  with 
published  Douglas-fir  height  growth  curves.  Because  the  published  curves  could  not 
adequately  describe  height  growth  patterns  of  these  trees,  new  height  growth  and  site 
index  curves  were  constructed. 

Four  estimating  equations  for  height  growth  and  three  for  site  index  were  fitted  by  use 
of  weighted  least  squares  regression,  and  site  index  and  height  were  constrained  to  be 
equal  at  index  age,  100  years.  The  equations  with  the  best  fit  to  the  data  across  the 
ranges  of  age  and  site  index  were  chosen  as  the  final  models.  The  resulting  curves  differ 
from  those  previously  published  for  Douglas-fir. 


Instructions  for  using  the  equations  are  presented. 
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Introduction  In  the  western  Cascade  Range  of  Oregon  and  Washington,  height  growth  and  site  index 

curves  are  available  for  Douglas-fir  (Pseudotsuga  menziesii  (Mirb.)  Franco)  on  mesic, 
midelevation  sites  (McArdle  and  others  1 961 ,  King  1 966),  and  on  upper  slopes  near  the 
cool  extreme  of  its  range  (Curtis  and  others  1974b).  No  height  growth  or  site  index 
curves  have  been  available,  however,  for  Douglas-fir  near  the  relatively  hot,  dry  extreme 
of  its  range  in  the  western  Cascades.  Height  growth  of  Douglas-fir  is  slower  initially  but 
is  more  prolonged  on  such  dry  sites  than  on  moister  sites  (Means  1 980),  and  as  a  result 
volume  growth  follows  a  similar  trend.  Thus,  available  site  index  curves  (for  example, 
McArdle  and  others  1961,  King  1966),  yield  tables  (for  example,  McArdle  and  others 
1 961 ,  Curtis  and  others  1 982),  and  height  growth  curves  in  stand  growth  simulators  (for 
example,  DFSIM  by  Curtis  and  others  1981)  are  not  applicable.  The  objective  of  this 
study  was  to  construct  height  growth  and  site  index  equations  for  Douglas-fir  on  hot,  dry 
sites  in  the  Willamette  National  Forest  in  Oregon. 

Curves  that  reach  well  beyond  200  years  are  needed  for  estimating  site  index  for  old 
stands  and  for  estimating  height  at  advanced  ages  to  evaluate  the  effects  of  manage- 
ment of  stands  for  old-growth  characteristics. 

Several  methods  have  been  used  in  constructing  Douglas-fir  height  growth  and  site 
index  curves.  McArdle  and  others  (1 961 )  built  height  growth  curves  from  many  one-time 
observations  of  height  and  age,  fitting  a  mean  guide  curve,  then  making  other  curves 
proportional  to  the  guide  curve.  This  technique,  however,  includes  several  sources  of 
significant  error  that  can  be  eliminated  by  using  stem  dissection  data  (Curtis  1 964).  King 
(1 966)  used  stem  dissection  data  to  construct  height  growth  curves  for  Douglas-fir  with 
height  as  the  dependent  variable  and  then  used  these  curves  to  estimate  site  index.  Site 
index  estimating  equations,  however,  should  be  fit  with  site  index  as  the  dependent 
variable,  and  height  growth  estimating  equations  should  be  fit  with  height  as  the 
dependent  variable  (Curtis  and  others  1974a).  Curtis  and  others  (1974b)  built  upper- 
slope  Douglas-fir  height  growth  and  site  index  curves  in  this  way.  Dahms  (1 975)  built  on 
this  basic  idea  by  first  fitting,  for  decadal  ages,  linear  regressions  to  estimate  site  index 
(SI  =  aO  +  a1  HT)  and  height  (HT  -  bO  +  b1  SI),  then  building  age  back  into  the 
models  by  defining  the  slopes  of  these  regressions  and  the  mean  height  curves  in  terms 
of  age;  Dahms  worked  with  lodgepole  pine  {Pinus  contorta  Dougl.  ex  Loud.).  Cochran 
(1979)  used  Dahms'  technique  to  build  curves  for  Douglas-fir  in  eastern  Oregon  and 
Washington.  Monserud  (1984)  found  Dahms'  technique  suitable  as  a  first  step  in 
deriving  a  site  index  estimating  equation,  but  inferior  to  a  logistic  equation  for  describing 
height  growth  of  Douglas-fir  in  the  northern  Rocky  Mountains.  As  a  final  step,  Monserud 
(1 984)  eliminated  highly  correlated  terms  from  the  site  index  equation  without  degrading 
model  fit.  The  equation  forms  used  for  Douglas-fir  by  the  above  investigators  were 
examined  in  this  study. 
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This  study  was  a  component  of  a  larger  study  in  the  Detroit,  Blue  River,  McKenzie, 
Oakridge,  and  Rigdon  Ranger  Districts  of  the  Willamette  National  Forest  (Means  1 980). 
Trees  were  cut  in  all  these  Districts  except  Detroit,  where  dry  coniferous  forest  is 
relatively  limited  in  extent,  so  suitable  sampling  sites  could  not  be  found  (fig.  1). 

Trees  were  cut  from  dry  coniferous  forests,  defined  as  climax  Douglas-fir  and 
incense-cedar  (Libocedrus  decurrens  Torr.)  and  lacking  significant  (0.1  percent  cover) 
western  hemlock  (Tsuga  heterophylla  (Raf.)  Sarg.)  in  any  size  class.  These  dry-site 
Douglas-fir  community  types  (Means  1980)  are  more  similar  to  the  Pseudotsuga/ 
Holodiscus  type  than  to  any  other  community  type  of  Dyrness  and  others  (1974). 

Tree  Selection  Several  criteria  were  used  for  selecting  trees.  Candidate  trees  were  dominants  or  strong 

codominants  with  straight  boles  and  no  sign  of  top  breakage.  Unhealthy  crowns,  flat  or 
deformed  tops,  or  several  conks  (or  one  large  conk)  were  causes  for  rejection. 
Candidates  were  increment  bored  at  breast  height.  Trees  showing  periods  of 
suppressed  radial  growth  or  significant  rot  were  not  used.  These  criteria  for  increment 
cores  resulted  in  rejection  of  70  to  90  percent  of  the  trees  more  than  200  years  old,  and 
30  to  70  percent  of  the  trees  100  to  200  years  old  that  met  the  exterior  criteria.  Effects 
of  past  fires  (scars  and  resultant  butt  rot),  top  breakage,  and  disease  increase  with  age. 
Suppressed  radial  growth  was  a  common  cause  for  rejection,  which  is  not  surprising 
because  canopy  trees  often  become  established  in  stands  with  mature  trees  (Means 
1982).  Many  plots  in  stands  with  canopy  trees  more  than  200  years  old  contained  no 
trees  suitable  for  dissection  (and  therefore  none  suitable  for  estimation  of  site  index). 
The  high  rate  of  rejection  precluded  assessment  of  the  potential  bias  in  site  index  curves 
caused  by  changes  in  the  tallest  tree  on  a  plot  through  time  (Dahms  1963). 

Stem  Dissection  The  stems  of  40  trees  on  thirty-six  0.1  -ha  plots  were  dissected ;  in  general,  the  techniques 

of  Herman  and  others  (1 975)  were  followed.  In  most  cases  cross-section  disks  were  cut 
at  stump  height,  at  breast  height  (1 37  cm  above  ground  level),  and  at  2-m  intervals  above 
breast  height.  Ground  level  was  defined  as  the  mean  of  the  ground  levels  on  the  uphill 
and  the  downhill  sides  of  the  tree.  In  1 977  whole  disks  were  taken  to  the  laboratory,  and 
narrow  (approximately  2x2  cm)  sections,  the  length  being  the  distance  from  the  edge  to 
the  pith  of  the  disk,  were  cut  with  a  table  saw.  Cutting  these  narrow  sections  precluded 
cross-checking  ring  counts  on  additional  radii.  Therefore,  in  1978  wider  sections  (10  to 
20  cm)  were  cut  from  disks  in  the  field  with  a  chain  saw;  and  whole  disks  were  usually 
collected  at  stump  and  breast  height. 

Data  Preparation  The  narrow  sections  from  trees  cut  in  1977  were  prepared  by  sanding.  The  wider  1978 

sections  were  prepared  with  a  belt  sander  if  the  wood  was  dry  or  with  a  "Surform"-  plane 

Ring  Counting  if  the  wood  was  moist.  Narrow  rings  (less  than  about  1  mm)  were  sanded  and  were 

counted  under  a  binocular  microscope  with  variable  magnification.  Maintaining  standing 
water  along  radii  improved  our  ability  to  distinguish  rings  but  was  seldom  done  because 
it  was  difficult  to  sustain  for  more  than  a  minute.  Usually  two  radii  were  counted  at  each 
disk  height. 

U  The  use  of  trade,  firm,  or  corporation  names  in  this  publication 
is  for  the  information  and  convenience  of  the  reader.  Such  use 
does  not  constitute  an  official  endorsement  or  approval  by  the 
U.S.  Department  of  Agriculture  of  any  product  or  service  to  the 
exclusion  of  others  that  may  be  suitable. 


Field  Methods 

Study  Area 
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Figure  1. — Locations  of  sample  trees  (•)  in  the  Willamette  National 
Forest  (RD  is  Ranger  District). 


I* 

Data  Cleaning  Height  versus  age  was  then  plotted  for  each  tree.  Anomalies  such  as  abrupt  changes  in  ^ 

and  Screening  height  growth  rate  were  often  corrected  by  recounts  or  counts  of  new  radii  on  whole 

sections.  A  program  written  in  FORTRAN  77,  STEM,  was  used  to  plot  height  over  age 
and  height  growth  rate  for  each  tree  and  to  calculate  (by  linear  interpolation)  heights  at 
decadal  ages  for  fitting  height  growth  and  site  index  models. 

Ten  trees  showed  suppressed  height  growth  early  in  life  or  had  isolated  periods  of 
relatively  slow  height  growth,  indicating  significant  top  breakage  or  suppression.  These 
were  excluded  from  further  analysis.  Two  trees  with  breast-height  age  less  than  100 
were  too  young  to  be  used  for  constructing  site  curves  with  an  index  age  of  100  and 
were  also  excluded.  After  these  trees  were  rejected,  only  one  0.1  -ha  plot  remained  with 
two  trees.  Height  growth  curves  for  these  trees  did  not  cross,  so  their  decadal  heights 
were  averaged  and  treated  as  one  tree.  The  height  growth  curves  of  two  trees  were 
truncated  to  eliminate  height  growth  anomalies  at  advanced  ages.  The  analyses  thus 
were  based  on  27  trees  and  505  observations  of:  height  above  breast  height,  site  index 
(height  above  breast  height  at  1 00  years  breast-height  age),  and  decadal  breast-height 
age.  These  trees  are  summarized  by  site  index  class,  breast-height  age,  and 
topographic  position  in  table  1,  and  by  dry-site  piant  community  in  table  2. 


3 


Table  1 — Distribution  of  the  27  Douglas-fir  trees  used  in  analyses,  by  breast-height 
age,  elevation,  slope,  and  aspect  for  ranges  of  site  index 




Range  of  site  index 
(meters  above  breast  height,   index  age  100) 


All  site 

Stratum  26-30      30.1-34    34.1-38     38.1-42    42.1-46    46.1-50  indices 


Number  of  trees 


Breast-height 
age  (yr): 
100-150 
151-200 
201-250 
251-300 


1  5 

1  2 

1  1 

1  2 


10 

3 
5 
9 


Years 


Mean  age 


Elevation  (m) 
500-599 
600-699 
700-799 
800-899 
900-999 


240  263  200  172  122 

------------  Number  of  trees 


112 


194 


1 

2 
4 
1 

2 

Meters 


3 
4 
12 
5 
3 


Mean 
elevation 


631  779  760  776  771 

------------  Number  of  trees 


899 


767 


Slope  (percent) 
0-20 
21-40 
41-60 
61-80 
81-100 


1 

1  1 

2  4  1 
4  1 

1 


Percent 


Mean  slope 


39 


Aspect  (degrees): 
91-180 
181-270 
2  70-360 

Total  trees 


62 


50  53  58 

■  -  -  Number  of  trees 


3 
6 
1 

10 


67 


54 


6 
19 

2 

27 


Table  2 — Distribution  of  the  27  Douglas-fir  trees  used  in  analyses,  by  plant 
community  for  ranges  of  site  index 


Range  of  site  index 
(meters  above  breast  height,  index  age  100) 


Plant  All  site 

community  1/  26-30      30.1-34    34.1-38    38.1-42    42.1-46    46.1-50  indices 


-----------  Number  of  trees  ---------- 

Psme/ Hodi/ grass : 

Asde  phase  1  1 

Cohe  phase  2                                11  4 

Psme /Hodi -Acci  1                                 12  12  7 

Psme/Beaq/Dispo  12                11  5 

Lide/Whmo                                                                       2  2 

Lide/ Chum                                   2                13  2  8 


1/  Plant  communities  are  described  by  Means  (1980).     The  species  codes  are: 
Acci  =  Acer  c  ircinatum,  Asde  =  Aspidot is  densa ,   Beaq  =  Berberis  aquifolium, 
Chum  =  Chimaphi la  umbel lata ,  Cohe  =  Col lomia  heterophylla,  Dispo  =  Disporum 
spp.,  Hodi  =  Holodiscus  discolor ,   Lide  =  Libocedrus  decurrens ,   Psme  =  Pseudo- 
tsuga  menzies  ii ,     Whmo  =  Whipplea  modesta. 


Analysis 

Examination  of 
Available  Curves 


Because  our  data  set  is  small  (27  trees)  it  was  important  to  determine  as  a  first  step 
whether  new  curves  were  really  needed  or  whether  existing  Douglas-fir  curves  would 
suffice.  Toward  this  end,  we  compared  our  height  growth  curves,  smoothed  using 
Heger's  (1 968)  data  summarization  method  (at  site  indices  of  28,  32,  36, 40,  and  44  m), 
with  the  height  growth  curves  published  by  Curtis  and  others  (1 974b),  Monserud  (1 984), 
King  (1966),  and  Cochran  (1979).  Site  indices  for  these  curves  were  chosen  so  they 
would  cross  our  curves  at  their  index  ages. 

The  height  growth  curves  of  King  (1966),  Cochran  (1979),  and  Monserud  (1984)  were 
consistently  biased  compared  with  ours;  they  had  higher  values  for  height  before  their 
index  age  (50  yr)  and  lower  values  past  this  age.  At  age  1 20,  discrepancies  were  2.5  to 
3.5  m  for  Cochran's,  4  to  6.5  m  for  King's,  and  3  to  7.5  m  for  Monserud's.  These 
discrepancies  consistently  increased  past  age  120,  especially  for  Cochran's  (1979) 
height  curves  which  declined  past  age  1 80-an  extrapolation  past  his  data  that  extended 
to  1 30  years. 

Curtis  and  others'  (1979b)  upper-slope  Douglas-fir  height  growth  curves  are  closer  to 
ours  in  a  narrow  range  of  moderate  site  indices  (32-36  m),  where  errors  are  all  less  than 
2  m  (within  the  range  of  our  data).  At  low  sites  indices,  however,  our  dry-site  trees  grow 
more  slowly  initially  and  more  rapidly  at  advanced  ages  than  do  upper-slope 
Douglas-firs.  At  high  site  indices,  conversely,  our  trees  grow  more  rapidly  initially,  and 
more  slowly  later  than  do  upper-slope  Douglas-firs.  Errors  range  from  -1 .5  to  2.2  m  for 
trees  less  than  1 00  years  old  and  from  -2.7  to  3.0  m  for  trees  older  than  1 00  years.  The 
dry-site  Douglas-fir  height  growth  curves  are  more  polymorphic  than  the  curves  for 
upper-slope  Douglas-fir. 
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Choice  of  Index  Age 


Because  height  growth  curves  derived  by  Heger's  (1968)  method  differed  from  those 
already  available,  our  site  index  curves  for  Douglas-fir  on  hot  dry  sites  would  also  be 
different;  therefore  we  decided  to  build  new  curves  even  though  our  data  set  was  small. 

An  index  age  of  100  years  at  breast  height  (1 .37  m)  was  chosen  because  this  will 
probably  be  near  rotation  age  in  managed  stands  of  this  type,  and  the  reliability  of  site 
index  estimates  should  be  better  than  it  would  be  if  a  younger  index  age  were  chosen. 
Site  index  estimates  are  less  accurate  when  differences  between  stand  age  and  index 
age  are  relatively  large,  and  many  dry  coniferous  forest  stands  are  1 00  to  250  years  old. 


Weighted  Regression 


All  height  growth  and  site  index  models  were  fitted  by  weighted  least  squares  regression 
since  variance  was  not  uniform  across  the  range  of  ages.  The  iterative  weighting 
technique  of  Monserud  (1984)  was  used,  in  which  the  weight  was  the  inverse  of  the 
variance  at  each  decadal  age. 


Constraining  the 
Models 


Height  growth  models  were  constrained  so  that  trees  "grown"  in  computer  models  will 
reach  a  height  equal  to  known  site  index  at  index  age.  For  all  (except  one)  height  growth 
models,  HT  =  f(SI,age),  height  (HT)  was  constrained  to  equal  site  index  (SI)  at  index 
age  by  subtracting  SI  =  f(SI, index  age)  from  the  height  model  before  fitting,  giving 
HT  =  SI  +  f(SI,age)  -  f(SI, index  age)  (the  model  of  Curtis  and  others  (1974b)  was 
constrained  using  the  technique  they  describe).  Site  index  models,  SI  =  f(HT,age),  were 
similarly  constrained  so  that  trees  measured  in  the  field  at  the  index  age  would  have 
current  height  as  their  site  index.  Errors  up  to  0.5  m  at  index  age  would  occur  if  equations 
were  not  constrained.  Constraining  typically  reduced  divergence  from  Heger's  (1968) 
curves  (see  below)  within  the  range  of  most  of  our  data  and  rarely  degraded  the  fit. 


Construction  of  Height 
Growth  Curves 


Heger's  curves. — The  first  step  in  construction  of  height  growth  curves  was  to  fit, 
individually  for  decadal  breast-height  ages,  linear  regressions  of  height  on  site  index 
(Heger  1968)  of  the  form 

HT  =  a  +  bSI;  (1) 


where  HT  =  total  height  - 1 .37  m  (height  above  breast  height),  SI  =  HT  at  1 00  yr,  and 
a,  b  =  regression  coefficients.  A  total  of  28  regressions  were  fit,  for  decadal  ages  1 0  to 
280.  Examination  of  the  data  and  residuals  showed  no  curvilinear  trends  at  any  age.  We 
constructed  height-over-age  curves  (fig.  2)  from  this  family  of  28  regressions  and 
compared  them  with  candidate  height  growth  curves  to  aid  in  model  selection.  Heger 
(1968)  used  such  curves  as  height  growth  curves,  but  we  derived  a  more  general 
equation  because  these  were  unrealistically  jagged  at  advanced  ages,  and  to  provide 
one  equation  for  use  in  computer  programs. 
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Figure  2. — Height  growth  curves  based  on  Heger's  (1968)  family  of 
regressions,  equation  (1) — solid  and  dashed  curves  (dashes 
indicate  points  beyond  the  data);  and  the  final  height  growth  model, 
equation  (6) — dotted  curves. 


Choosing  the  height  growth  model. — Four  height  growth  models  were  examined.  The 
model  used  by  Curtis  and  others  (1974b),  rewritten  to  express  height  as  a  function  of 
site  and  age,  is 


HT  =  Sl/(a  +  bAn  +  cAn/SI); 


(2) 


where  a,  b,  c,  n  =  regression  coefficients  and  A  =  breast-height  age.  The  model  King 
(1966)  used  was  examined, 


HT  =  A2/(a  +  bA  +  cA2); 
where  a,  b,  c  =  functions  of  SI.  The  equation  used  by  Beck  (1971)  is 

HT  =  a  (1  -  exp(bA))c; 


(3) 


(4) 


where  a,  b,  c  =  functions  of  SI  and  exp(bA)  is  the  base  of  the  natural  logarithm  to  the 
power  bA.  The  model  form  used  by  Monserud  (1984)  was  also  fit, 


HT  =  a  Slb/(1  +  exp(c  +  d  InA  +  e  InSI); 


(5) 


where  a,  b,  c,  d,  e  =  regression  coefficients. 
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Each  model  was  first  fit  to  the  smoothed  height  over  age  curves  derived  by  Heger's 
(1968)  method,  equation  (1).  At  all  but  the  lowest  site  index  levels,  model  (5) 
underestimated  height  at  low  ages  and  overestimated  it  at  high  ages  so  it  was  not 
examined  further.  The  remaining  models  were  fit  to  the  decadal  heights  and  ages,  and 
they  gave  nearly  identical  fits — standard  errors  of  the  estimates  (SEE's)  ranged  from 
1.73  to  1.76  m.  Examination  of  residuals  showed  no  problems  with  bias.  Equation  (2) 
was  chosen  because  it  had  a  simpler  form;  the  coefficients  of  equations  (3)  and  (4)  were 
functions  of  site  index  with  two  or  three  terms  each.  Its  final  form  after  it  was  constrained 
(Curtis  and  others  1974b)  was, 

HT  =  Sl/(d  +  b  An  +  e/Sl  +  c  A7SI);  (6) 

where  HT,  SI,  and  A  are  defined  as  for  equations  (1 )  and  (2);  d  =  0.90031 ;  b  =  36.932; 
e  =  -1 1 .9659;  c  =  4433.1 ;  and  n  =  -1 .28438;  with  an  SEE  of  1 .73  m.  Height  growth 
curves  based  on  equation  (6)  and  on  Heger's  regressions,  equation  (1),  are  shown  in 
figure  2. 

Construction  of  Site  Heger's  curves. — The  first  step  in  construction  of  site  index  curves,  similar  to  that  for 

Index  Curves  height  growth  curves,  was  to  use  Heger's  (1968)  method  to  fit,  individually  for  decadal 

ages,  linear  regressions  of  site  index  on  height  of  the  form, 

SI  =  a  +  b  HT;  (7) 

where  a,  b  =  regression  coefficients.  Examination  of  the  data  and  residuals  showed  no 
curvilinear  trends  at  any  age.  Inverse  estimates  of  height  for  several  levels  of  site  index 
were  calculated  and  used  in  plotting  these  relationships  in  the  traditional  height  over  age 
format.  The  linear  regressions,  equation  (7),  and  the  height  over  age  curves  were 
compared  graphically  with  candidate  site  index  estimating  equations  to  aid  in  model 
selection. 

Constructing  and  choosing  the  site  index  model. — Three  site  index  estimating 
equations  were  examined.  The  first  model  was  a  generalized  form  of  equation  (7)  (Curtis 
and  others  1974b), 

SI  =  a  +  b  HT;  (8) 

where  a,  b  =  functions  of  breast-height  age.  Stepwise  regression  was  used  to  select 
from  many  independent  variables  for  describing  a  and  b.  Monserud  (1 984)  derived  a  site 
index  model  by  first  constructing  an  equation  by  Dahms'  (1 975)  method,  then  removing 
highly  correlated  terms  and  recalculating  least  squares  parameter  estimates.  His  final 
model  was  the  second  model  examined  here: 

SI  =  a  +  b  (InA)2  +  c  A  InA  +  d  HT  +  e  HT  A;  (9) 

where  a,  b,  c,  d,  e  =  regression  coefficients.  The  third  model,  constructed  by  Dahms' 
(1 975)  method,  involves  finding  a  function  of  breast-height  age  (A)  for  b  in  equation  (8), 

b  =  fb(A),  (10) 

and  fitting  an  equation  to  mean  decadal  height  (MHT), 

MHT  =  fmht(A).  (11) 
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When  fitting  equation  (1 1),  we  wished  to  maintain  MHT  at  the  same  mean  site  index 
(MSI  =  37.31 2  m)  for  all  decadal  ages.  Therefore,  when  trees  dropped  out  at  advanced 
ages,  we  estimated  MHT  from  Heger's  regressions,  equation  (7),  at  SI  =  MSI.  Then 
MHT  was  substituted  for  HT,  and  MSI  (MSI  can  be  considered  constant  with  age  since 
older  MHT's  were  adjusted  to  it)  was  substituted  for  SI  in  equation  (8),  which  was  solved 
to  provide  an  equation  for  a: 

a  =  MSI  -  b  MHT.  (12) 

Substituting  equation  (11)  for  MHT  in  equation  (12),  then  substituting  equation  (10)  for 
b  and  the  modified  equation  (12)  for  a  in  equation  (8),  and  simplifying,  gave  a  smoothed 
site  index  estimating  equation  of  the  type  used  by  Dahms  (1975): 

SI  =  MSI  +  fb(A)  (HT  -  fmht(A)).  (13) 

This  equation  was  constrained  to  pass  through  HT  =  SI  at  index  age,  which  eliminated 
the  constant  MSI;  and  new  parameters  were  estimated  by  least  squares  to  improve  its  fit. 

Height  over  age  curves  for  each  site  index  model  were  compared  with  Heger's  curves 
from  equation  (7),  and  plots  of  the  "residuals"  (difference)  between  each  model  and 
Heger's  curves  were  examined.  Plots  of  site  (predicted,  actual,  and  Heger's)  over  height 
and  their  residuals  for  alternate  decadal  ages  were  also  examined.  Problems  with  bias 
became  readily  apparent,  which  facilitated  rejection  of  some  models  and  simplified 
selection  of  equation  (13)  as  the  final  model. 

Attempts  to  reduce  the  number  of  coefficients  of  equation  (13),  as  Monserud  (1 984)  did, 
markedly  reduced  the  fit,  so  all  fitted  coefficients  (b-g)  were  retained.  After  equation  (13) 
was  constrained  and  simplified  it  became 

SI  =  (b  +  c  A  +  d  A"0  5)  (HT  -  1  (e  +  f  A9))  +  k  HT  +  j;  (14) 

where  SI,  HT,  and  A  are  defined  as  for  equations  (1)  and  (2),  b  =  0.94622, 
c  =  -0.00175322,  d  =  3.6660,  e  =  0.0168601,  f  =  2.7676,  g  =  -1.19554, 
k  =  -0.137502,  and  j  =  40.470.  Site  index  curves  based  on  this  equation  are 
presented  in  figure  3  and  are  compared  with  those  based  on  Heger's  regressions, 
equation  (7),  in  the  more  familiar  inverted  form  in  figure  4. 


The  final  height  growth  and  site  index  curves  are  compared  in  figure  5. 
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Figure  3. — Final  site  index  estimating  curves,  equation  (14), 
presented  with  the  dependent  variable  (site  index)  on  the  vertical 
axis. 
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Figure  4. — Site  index  curves  based  on  Heger  s  (1968)  family  of 
regressions,  equation  (7) — solid  and  dashed  curves  (dashes 
indicate  points  beyond  the  data),  and  the  final  site  index  model, 
equation  (14) — dotted  curves. 
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Figure  5. — Comparison  of  final  height  growth  (solid  curves)  and 
site  index  curves  (dashed  curves)- 


Discussion  The  final  height  growth  curves,  equation  (6),  are  compared  with  other  height  growth 

curves  for  Douglas-fir  in  figure  6.  The  curves  of  McArdle  and  others  (1961)  for  western 
Comparison  With  Other      Washington  and  Oregon,  King  (1966)  for  Washington,  Cochran  (1979)  for  eastern 
Height  Growth  Curves        Washington  and  Oregon,  and  Monserud  (1984)  for  the  northern  Rocky  Mountains  all 

show  more  rapid  height  growth  at  early  ages  and  slower  height  growth  later  than  do  the 
curves  for  dry-site  Douglas-fir.  One  important  implication  is  that  volume  growth  in 
dry-site  Douglas-fir  stands  is  initially  slower  but  is  more  prolonged  and  maximum  mean 
annual  increment  occurs  later  than  in  the  other  forest  types  studied  (Means  1980). 

Height  growth  of  upper-slope  Douglas-fir  in  the  Washington  and  Oregon  Cascade 
Range  (Curtis  and  others  1 974b)  is  similar  to  that  on  dry  sites  at  midsite  indices  (fig.  6a). 
At  higher  site  indices,  however,  height  growth  on  dry  sites  is  initially  faster  but  is  slower 
later  in  life,  whereas  at  lower  site  indices  it  is  initially  slower  but  is  faster  later  in  life. 
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Figure  6. — Comparison  of  our  height  growth  curves,  equation  (6), 
with  height  growth  curves  previously  published  for  Dougias-fir.  Site 
indices  for  these  curves  were  chosen  so  they  would  cross  our  curves 
at  their  index  ages:  A.  companson  of  our  curves  with  those  of  Curtis 
and  others  (1974b)  ano  McArdie  ana  others  (1961 };  B.  comparison 
of  our  curves  with  those  of  Cochran  (1979).  King  (1966).  and 
Monserud  (1984). 


Reliability  The  reliability  of  all  predicting  equations  is  strongly  influenced  by  the  unexplained 

variation  in  the  data  used  to  build  the  equations.  The  variation  the  final  models  cannot 
explain  is  expressed  as  the  standard  error  of  the  estimate  (SEE)  in  figures  7  and  8. 
Estimates  are  less  reliable  when  the  SEE  is  greater. 


0  50  100  150  200  250  300 


Breast  height  age  (years) 

Figure  7. — Standard  error  of  the  estimate  (uncorrected  for  the 
degrees  of  freedom  in  the  models)  for  height  growth  models  by 
decadal  age.  The  solid  line  is  for  the  final  height  growth  model, 
equation  (6),  and  the  dashed  line  is  for  the  regressions  produced  by 
Heger's  (1968)  method,  equation  (1). 
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Figure  8. — Standard  error  of  the  estimate  (uncorrected  for  the 
degrees  of  freedom  in  the  models)  for  site  index  models  by  decadal 
age.  The  solid  line  is  for  the  final  site  index  model,  equation  (14),  and 
the  dashed  line  is  for  the  regressions  produced  by  Heger's  (1968) 
method,  equation  (7). 
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Reliability  is  reduced  by  use  of  the  equations  beyond  the  range  of  the  data  used  to 
construct  them.  This  range  is  indicated  in  figures  2  and  4.  Note  that  no  trees  older  than 
139  years  at  breast  height  were  found  at  high  site  indices  (table  1).  This  is  reflected  in 
figures  2  and  4  by  the  jagged  appearance  of  Heger's  curves  and  their  divergence  from 
the  final  model  curves  at  these  ages  and  sites.  As  a  result,  estimates  at  advanced  ages 
and  high  sites  are  relatively  unreliable  extrapolations. 

The  equations  presented  here  have  SEE's  similar  to  those  of  Monserud  (1984)  for 
Douglas-fir  in  the  northern  Rocky  Mountains,  the  only  other  Douglas-fir  curves  for  which 
this  statistic  was  found.  Our  SEE's  for  height  growth  equation  (6)  and  site  index  equation 
(14)  are  1.73  and  2.07  m,  respectively;  Monserud's  "tree-about-the-model"  errors  are 
1 .55  and  1 .89  m,  respectively.  When  data  for  ages  older  than  200  years  are  excluded 
(as  Monserud  did),  our  SEE's  (1 .53  and  2.06  m)  for  dry-site  Douglas-fir  are  even  closer 
to  his.  We,  however,  included  these  data  to  provide  a  better  basis  for  predictions  at 
advanced  ages,  and  because  the  data  for  these  ages  received  low  weights  in  the 
weighted  regression  procedure. 

We  considered  Heger's  linear  equations  (1 )  and  (7)  to  provide  the  best  possible  fit,  since 
there  was  no  curvilinear  relationship  between  site  and  height  at  any  age.  The 
uncorrected  sums  of  squares  for  error  (SSE)  of  the  final  models  were  similar  to  those 
for  Heger's  model  (figs.  7  and  8)  at  all  but  the  oldest  ages.  Better  final  models,  therefore, 
could  offer  only  slight  improvement  unless  new  variables  (for  example,  soils,  nutrients) 
were  considered.  We  sampled  a  narrow  range  of  plant  communities:  the  range  in  site 
index  (28.57  to  48.00  m),  only  52  percent  of  the  mean  site  index  (37.31  m),  is  lower  than 
that  of  other  studies.  Thus,  environmental  variables  are  unlikely  to  improve  model  fit. 

Overfitting  the  data,  constructing  models  that  describe  meaningless  irregularities  in  the 
data,  can  be  a  problem  when  models  are  fit  to  small  data  sets  such  as  this  one  of  27 
trees.  If  this  had  occurred,  it  is  reasonable  to  conclude  that  the  SEE's  would  have  been 
artificially  low.  The  previous  comparison  with  Monserud's  (1984)  SEE's,  however, 
indicates  that  this  probably  did  not  occur. 

Occasionally,  straight  boles,  well-formed  tops,  and  increment  cores  indicated  no 
problems  in  height  growth,  but  after  the  stems  were  dissected,  height  growth  curves  or 
portions  thereof  showing  height  growth  suppression  were  excluded  from  the  analysis. 
This  will  not  be  possible,  however,  in  routine  applications  of  these  equations;  in 
particular,  estimates  of  site  index  for  such  trees  will  be  low  and  estimates  of  height 
growth  (when  site  index  is  known  from  other  trees)  will  be  high.  Applications  of  the 
equations  in  young  stands  will  not  have  this  bias  because  all  the  rejected  trees  were  1 1 7 
years  at  breast  height  or  older.  This  problem  is  common  (though  seldom  discussed)  to 
all  studies  in  which  data  are  rejected  based  on  examination  of  height  growth  curves  after 
stem  dissection.  Examination  of  the  height  growth  curves  for  trees  rejected  in  this  study 
indicates  errors  in  site  index  and  height  will  usually  be  less  than  2-3  m. 

This  source  of  error  will  probably  be  more  important  in  dry-site  Douglas-fir  communities, 
where  potential  site  index  trees  often  regenerate  after  wildfire  in  stands  with  surviving 
trees  (Means  1 982),  and  will  be  less  important  in  community  types  that  are  usually  even 
aged.  Therefore  site  index  may  be  less  suitable  as  a  productivity  index,  and  the 
availability  of  growth  models  suited  to  uneven-aged,  dry-site  stands  is  probably  more 
critical. 


How  To  Use  the  The  height  growth  and  site  index  curves  differ  in  form  (fig.  5)  and  recommended  use. 

Equations  The  height  growth  equation  was  fit  with  height  as  the  dependent  variable.  It  describes 

the  height  growth  of  trees  that  attain  a  specified  site  index  at  index  age  and  is 
appropriate  for  prediction  of  height  and  height  growth  for  stands  of  given  site  index,  as 
in  construction  of  yield  tables  and  in  stand  simulators.  The  site  index  equation  was  fit 
with  site  index  as  the  dependent  variable.  It  is  a  more  efficient  estimator  of  site  index 
than  is  the  height  growth  equation  (Curtis  and  others  1974a)  and  should  be  used  for 
estimating  site  index  when  age  and  height  are  known. 

Stand  Selection  These  equations  should  be  applied  to  stands  in  the  Willamette  National  Forest  in  the 

Ranger  Districts  sampled  (fig.  1 ).  These  equations  may  be  more  appropriate  than  other 
available  equations  for  use  in  other  Ranger  Districts  or  in  adjacent  areas  (for  example, 
southwestern  Oregon);  this,  however,  is  not  certain.  The  error  associated  with  this  sort 
of  extrapolation  probably  increases  with  distance  from  the  areas  sampled  (fig.  1)  and 
cannot  be  estimated  from  data  at  hand. 

Stands  should  be  in  dry-site  Douglas-fir,  defined  as  Douglas-fir  or  incense-cedar  climax 
and  should  lack  a  significant  cover  (0.1  percent)  of  western  hemlock  in  any  size  class. 
The  above  comparisons  with  other  curves  indicate  that  the  dry-site  Douglas-fir  curves 
are  not  appropriate  for  other  communities. 

Tree  Selection  Site  index  estimates  should  be  based  on  one  or  two  (when  possible)  healthy  dominants 

or  strong  codominants  on  each  0.1 -ha  (0.25-acre)  plot.  Boles  should  be  straight  and 
there  must  be  no  signs  of  past  top  damage.  Candidate  trees  must  be  increment  bored 
at  breast  height,  and  trees  showing  suppressed  radial  growth  must  be  rejected.  The 
desired  precision  of  the  site  index  estimate  will  determine  the  number  of  plots  to  be  taken 
in  a  stand. 

Top  breakage  and  periods  of  suppression  are  likely  to  be  common  problems  on  dry  sites 
where  understory  trees  can  eventually  reach  the  canopy.  When  two  acceptable  trees 
occur  on  a  plot,  their  site  index  estimates  should  be  averaged.  Acceptable  site  index 
trees  are  rare  in  stands  with  many  trees  older  than  about  200  years,  so  estimating  the 
site  index  for  these  stands  may  be  difficult. 

Equations  for  Total  The  following  equations  are  in  a  form  probably  most  useful  to  most  people.  The  height 

Height  in  Feet  growth  equation,  equation  (6),  in  units  of  feet  is, 

HT  =  4.5  +  (SI-4.5)  /  (d  +  bAn  +  e/(SI-4.5)  +  cAn  (SI-4.5));  (15) 

where  HT  =  total  height  in  feet,  A  =  age  at  breast  height  (4.5  ft),  SI  =  total  height  in 
feet  at  1 00  yr,  d  =  0.90031 ,  b  =  36.932,  e  =  -39.258,  c  =  1 4544.39,  and  n  =  -1 .28438. 

The  site  index  estimating  equation,  equation  (14),  in  units  of  feet  is; 

SI  =  4.5  +  3.28084((b  +  cA  +  d/A°  5)(0.3048(HT-4.5)  -  1/(e  +  fA9)) 

+  k(HT-4.5)  +  j);  (16) 

where  SI,  HT,  and  A  are  defined  as  for  equation  (15),  b  =  0.94622,  c  =  -0.00175322, 
d  =  3.6660,e  =  0.0168601,  f  =  2.7676,  g  =  -1.19554,  k  =  -0.04191 1 ,  and  j  =  40.470. 
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Equations  and  curves  are  presented  for  estimating  height  and  site  index 
of  Douglas-fir  {Pseudotsuga  menziesii  (Mirb.)  Franco)  on  hot,  dry  sites  in 
the  Willamette  National  Forest  in  western  Oregon.  The  equations  are 
based  on  the  dissected  stems  of  27  trees.  The  curves  differ  from  those 
previously  published  for  Douglas-fir.  Instructions  are  presented  for  their 
use. 

Keywords:  Increment  (height),  site  index,  Douglas-fir,  Oregon  (Willamette 
National  Forest). 
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